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A massive white dwarf member of the Coma Berenices Open Cluster 
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ABSTRACT 



We report the identification, from a photometric, astrometric and spectroscopic study, 
of a massive white dwarf member of the nearby, approximately solar metalicity. Coma 
Berenices open star cluster (Melotte 111). We find the optical to near-IR energy distribu- 
tion of WD 12 16+260 to be entirely consistent with that of an isolated DA and determine 
the effective temperature and surface gravity of this object to be rcff=15739l}ggK and 



log 5=8.46: 



-0.03 
0.02- 



We set tight limits on the mass of a putative cool companion, M^O.O36M0 



(spatially unresolved) and M<;O.O34M0, (spatially resolved and a^2500AU). Based on the 
predictions of CO core, thick-H layer evolutionary models we determine the mass and cool- 
ing time of WD1216+260 to be Mwd=0.90 ± O.O4M0 and Tcooi=363i:^'jMyrs respectively. 
For an adopted cluster age of T=500±100Myrs we infer the mass of its progenitor star to be 
Minit=4. 77^0 97^0. We briefly discuss this result in the context of the form of the stellar 
initial mass-final mass relation. 

Key words: stars: white dwarfs; galaxy: open clusters and associations: Melotte 111, Coma 
Berenices open cluster 



1 INTRODUCTION 

Galactic open star clusters have for decades been regarded as key 
entities with which to address important issues in stellar and galac- 
tic astrophysics e.g the form of the stellar initial mass function, the 
evolution of stellar angular momentum and the form of the initial 
mass-final mass relation. This is because fundamental stellar pa- 
rameters such as age and composition can be more stringently con- 
strained through studying a co-eval ensemble of objects at essen- 
tially the same distance (e.g. Meynet, Mermilliod & Maeder 1993). 

Melotte 111 (Coma Berenices cluster, RA=12''22'", 
DEC=+26°, J2000.0) is the second closest open star cluster to 
the Sun. It appears to be marginally younger than the Hyades 
or Praesepe, with most recent age estimates lying in the range 
r=400-600Myrs (e.g. r^450Myrs, Bounatiro & Arimoto 1992; 
r^500Myrs, Odenkirchen et al. 2001; r?5i600Myrs, Kharchenko 
et al. 2005). Cayrel de Strobel (1990) has determined [Fe/H]=- 
0.065±0.021 from a high resolution spectroscopic study of 
a sample of eight F-,G- and K- type members. A similar but 
independent investigation of fourteen F and G dwarf associates 
has concluded that [Fe/H]=-0.052±0.047 (Friel & Boesgaard 
1992). More recenfly, Gebran, Monier & Richardo (2008), have 
determined [Fe/H]=0.07±0.09 through the high resolution spec- 
troscopic study of eleven F dwarf members. Thus the cluster 
metalicity is relatively close to the solar value. 
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In the last decade there have been several attempts to deter- 
mine the distance of Melotte 111. For example, using Hipparcos 
parallax measurements. Van Leeuwen (1999) and Robichon et al. 
(1999) have estimated d=89.9±2.1pc and d=87.0±1.6pc, respec- 
tively. However, the reliability of both of these measurements has 
been challenged by Pinsonneault et al. (1998) and Makarov (2003). 
The former matched a theoretical isochrone to the cluster sequence 
in the B-V colour to obtain d=80.9±1.5pc or d=84.1pc assum- 
ing [Fe/H]=-0.07±0.02 or solar metalicity respectively. The latter 
re-assessed the Hipparcos intermediate astrometry data using a 
method which reduced the impact of errors in the along-scan at- 
titude parameters to determine d=80.6±1.6pc. However, from his 
new reduction of the raw Hipparcos data, which included an 
improved treatment for the along-scan attitude issues and which 
is claimed to be superior to the Makarov (2003) approach, van 
Leeuwen (2007) has determined d=86.7lj!, gpc. This is consistent 
with the Geneva photometry based estimate of Nicolet (1981), 
d=85.4±4.9pc and is closer to the Pinsonneault et al. (1998) deter- 
mination, particularly if the cluster metalicity is nearer solar than 
originally assumed by these authors. As a consequence of this close 
proximity and the high galactic latitude of Melotte 111, levels of 
foreground reddening are very low e.g. E(B — 'l/)=0.006±0.013, 
Nicolet (1981); E{B - 1^)^^0.0032, Taylor (2006). 

The mean heliocentric radial velocity of cluster members is 
close to rv~Okms~^. For example, Abt & Willmarth (1999) mea- 
sured rv=-0.4±0.4 kms^^ from sixteen stars. Kharchenko et al. 
(2005) determined rv=-1.17±0.99 kms"^ by cross-referencing the 
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Figure 1. Left: Jtgr colour-colour diagram of stellar objects lying within the projected core radius of Melotte 111. Theoretical colour-colour tracks for DA 
white dwarfs with cooling times Si6()()Myrs, from Holberg & Bergeron (2006), are overplotted (thin grey line, O.7M0, medium grey line, O.9M0, thick grey 
line, I.IM0). Probable white dwarfs are highlighted (triangles). Right: g,u-i colour-magnitude diagram of the stellar objects lying within the projected core 
radius of Melotte 111. Theoretical white dwarf tracks are overplotted (as previously denoted). The locations in colour-magnitude space of the four probable 
white dwarfs (triangles) indicate that they are potentially members of Melotte 111. 



ASCC-2.5 catalogue (Kharchenko et al. 2001) with the General 
Catalogue of Radial Velocities (Barbier-Bros.sat & Figon 2000). 
From CORAVEL spectroscopic observations of the G7 III -1- 
A2.5rV giant binary member Trumpler 91 (Trumpler 1938), Mer- 
milliod, Mayor & Udry (2008) cstiinatc Tv =-l-1.06±l. Skins" ^ De- 
spite being nearby, the tangential motion of members is rather 
small. From Hipparcos data Robichon et al. (1999) determine 
Hacos 5=-11.38±0.23mas yr'S /i5=-9.05±0.12 mas yr'^ This 
is corroborated by the location of the clump of stars brighter than 
J=ll seen in a UCAC2 proper motion vector point diagram for 
the region of sky centered on the cluster (Figure 1 of Casewell et 
al. 2006). Kraus & Hillenbrand (2007) have recently performed 
a detailed study of the structure and the mass of Melotte 111 
utilising both astrometry and photometry. They have concluded 
that the projected tidal radius is 4.3±0.2°, which translates to 
6.4±0.3pc at a distance of d=85pc, and that the total cluster mass 
isMtot=112±16Mo. 

To date, there have been very few efforts to identify and study 
the white dwarf population of Melotte 111. This is presimiably be- 
cause the cluster is relatively sparse and its members do not have 
a distinctive proper motion. From a UBV photographic survey of 
21.2 sq. degrees of sky centered on the cluster down to mpg«17, 
Stephenson (1960) identified five candidate white dwarf members, 
but none of these appears to have yet been subject to detailed 
follow-up study. Brosch et al. (1998) used the FAUST ultraviolet 
imaging telescope to observe a 67.5 sq. degree region of sky cen- 
tered on the Coma galaxy cluster (RA=12'*58'", DEC=-i-28°05"', 
J2000.0) and claims to have unearthed ~5-10 candidate white 
dwarf members of Melotte 111. However, these objects all lie well 
beyond the tidal radius of the cluster and thus are clearly not grav- 
itationally bound members. 

Nevertheless, based on scaling the number of known white 
dwarf members of the slightly older, more massive Hyades 
(Mtot=3OO-46OM0, Pels et al. 1975, Reid 1992; NwdwIO, von 
Hippel 1998) and Praesepe (Mtot=500-600MQ, Adams et al. 2002, 



Kraus &. Hillenbrand 2007; Nwd~10, Casewell et al. 2008) clus- 
ters, we should expect a photometric survey of Melotte 1 1 1 to un- 
earth ~2 degenerate members. The detailed investigation of the 
white dwarf members of open star clusters (e.g. Sweeney 1976, 
Weidemann 1977, Wcidcmann 2000, Williams ct al. 2004, Kalirai 
et al. 2005, Dobbie et al. 2006a) is arguably the best observational 
approach to ascertaining the form of the stellar initial mass-final 
mass relation (IFMR). The difference between the age of an open 
star cluster and the cooling time of a white dwarf member provides 
an estimate of the lifetime of the progenitor star. The mass of the 
progenitor star can then be estimated by comparing this lifetime 
to the predictions of stellar evolutionary calculations. Robust con- 
straints on the form of the IFMR are very useful to several threads 
of astrophysical research. For example, the relation is an impor- 
tant component of models of the chemical evolution of the Galaxy 
as it provides an estimate of the amount of gas, eirriched with C, 
N and other metals, a low or intermediate mass star returns to the 
ISM (e.g. Carigi, Colin & Peimbert 1999). Understanding the form 
of the IFMR is crucial to deciphering information locked up in the 
white dwarf luminosity functions of stellar populations (e.g. Jeffery 
et al. 2007, Oswalt et al. 1996). Furthermore, the shape of the upper 
end of the IFMR is of particular interest as it places limits on the 
minimum mass of star that will experience a core-collapse super- 
nova explosion. With a tight handle on this mass, for example, the 
observed diffuse SNe neutrino background can serve as an empir- 
ical normalisation check on estimates of the star formation history 
of the universe (e.g. Hopkins et al. 2006). 

Driven by a desire to improve current observational con- 
straints on the form of the IFMR and the availability of high quality 
multiband photometry from the Sloan Digitial Sky Survey (SDSS) 
for Melotte 1 1 1, we have undertaken a search of the central regions 
of the cluster for white dwarf members. In the next section we out- 
line our photometric and astrometric survey of the cluster center. 
We then describe the acquisition and analysis of our spectroscopic 
follow-up data and the use of this to assess the membership status 
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Table 1. A summary of the astrometric and photometric properties of the four candidate white dwarf members of Melotte 111 unearthed by our 
SDSS based imaging survey of the central cluster region. 



Object 


Stephenson 


u 


9 


r 


i 


HaCOS S 






ID No. 




SDSS 




(mas yr~^) 


(mas yr~^) 


SDSS J121635.01+262934.7 


1 


14.88±0.03 


15.26±0.02 


15.78±0.02 


16.16±0.01 


-29.0±10.3 


-13.2±9.2 


SDSS J121845.69+264831.7 


2 


16.32±0.01 


16.47±0.01 


16.86±0.02 


17.20±0.03 


61.2±58.8 


-55.4±57.0 


SDSS J121856.17+254557.1 


3 


16.69±0.01 


16.52±0.02 


16.75±0.01 


17.()3±0.02 


-9.7±7.9 


-8.2±8.2 


SDSS J122420.51+264738.9 




16.95±0.03 


17.03±0.02 


17.36±0.02 


17.59±0.03 


2.5±6.2 


-10.0±6.1 



of candidates. Finally, we examine a likely white dwarf member in 
the context of the IFMR. 



2 A SEARCH FOR POTENTIAL WHITE DWARF 
MEMBERS OF MELOTTE 111 

2.1 Optical photometry 

Melotte 111 lies entirely within the footprint of the SDSS data re- 
lease 6 (DR6), so for this region of sky we have at otir disposal opti- 
cal photometry in five bands, ugriz, with a precision of up to 2-3% 
(Adelman-McCarthy et al. 2008). Harris et al. (2003) have demon- 
strated that single white dwarfs with Teff^l2000K are clearly sep- 
arated from main sequence objects and quasars in a ugr colour- 
colour diagram. Therefore, as a first step, we have extracted from 
the SDSS DR6 Data Archive Server, clean u,g and r band photom- 
etry for objects classified as stellar lying within 1.5° of the cluster 
center. This corresponds roughly to the projected e-folding radius 
for the surface density of O.S-l.OM© members (Kraus & Hillen- 
brand 2007). We have imposed a brightness criterion of 5^17.5 
since given the age and the distance of the cluster we would not ex- 
pect to find white dwarf members below this magnitude limit (e.g. 
Bergeron et al. 1995). The ugr colour-colour diagram for these data 
is shown in Figure 1, left. Model colours (Bergeron et al. 1995, Hol- 
berg & Bergeron 2006) for 0.7, 0.9 and I.IMq DA white dwarfs, 
with rcooi=^600Myrs, are overplotted (thin grey, medium grey and 
thick grey line respectively). Inspection of this plot reveals that four 
objects clearly stand out as candidate white dwarfs (filled trian- 
gles). Nevertheless, from this analysis alone we caimot determine 
if these are cluster members since field white dwarfs and some hot 
subdwarfs can also reside in this part of colour-colour space. In an 
attempt to exclude field contaminants we have examined their loca- 
tion in the u-i,g colour-magnitude diagram (Figure 1, right). Once 
again we have overplotted synthetic photometric tracks (Bergeron 
et al. 1995, Holberg & Bergeron 2006) for 0.7, 0.9 and I.IMq DA 
white dwarfs with TcooI ^ 600Myrs, here placed at d=85pc (denoted 
as in the colour-colour plot). However, given the effects of the fi- 
nite depth of Melotte 111 (rtidai=6.4±0.3pc) and the proximity of 
the four objects to these model tracks, at this stage we are led to re- 
tain all as candidate white dwarf cluster members. We note that our 
three brightest objects were previously identified as potential mem- 
bers of Melotte 1 1 1 by Stephenson (1960). The coordinates and the 
ugri photometry for all four candidates can be found in Table 1. 



2.2 Proper motions 

As a further aid to discriminating between cluster members and 
the field population, we have extracted proper motions for the four 
candidate white dwarfs from the SuperCOSMOS Science Archive 
(Table 1). These measurements, which are plotted in Figure 2, 




COS 6 (mas yr ^) 

Figure 2. A proper motion vector point diagram, based on SuperCOSMOS 
data obtained from the Survey Science Archive, for the photometrically se- 
lected candidate white dwarf members of Melotte 111. The mean proper 
motion of cluster members, as determined by Robichon et al. (1999) us- 
ing Hipparcos data, is highlighted (' x ')• As the proper motion of SDSS 
J122420.51+264738.9 lies ^2.2ct from the cluster mean we conclude that 
it is likely to be a field star. The proper motions of the remaining three po- 
tential white dwarfs can be considered at least marginally consistent with 
Melotte 111. 



have been determined with respect to background galaxies and, 
like the Hipparcos astrometry (e.g. Robichon et al. 1999), are ef- 
fectively absolute (see Hambly et al. 2001 for further details). Al- 
though the small proper motion of the cluster dictates that these 
SuperCOSMOS data are not of sufficient accuracy to compellingly 
demonstrate that a particular candidate is an associate of Melotte 
111, an object must have astrometric properties at least consis- 
tent with the cluster to be considered a member. Indeed, since 
SDSS I122420.51-I-264738.9 has a proper motion lying ^2.2cr 
from the Hipparcos determination of the cluster mean (' x ' in Fig- 
ure 2) we conclude that it is likely to be a field star and conse- 
quently do not consider it further in our analysis. While SDSS 
J121856.17-K254557.1 is the only candidate with a proper motion 
within 1(7 of that of Melotte 111, the measurements for SDSS 
1121635.01-^262934.7 and SDSS J121845.69-H264831.7 are at least 
marginally consistent with that of the cluster and do not rule out 
the possibility that either is a member. However, a recent spec- 
troscopic analysis of SDSS J121635.01-I-262934.7 (PG1214-I-268) 
has confirmed it to be a DA white dwarf with M~O.6M0 and 
Teff«65000K, residing at d~315pc (Liebert, Bergeron & Holberg 
2005). Thus we can say with confidence that it is not a member of 
Melotte 111 and exclude it from further discussion. Clearly, spec- 
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troscopic data is required to examine more thoroughly the nature 
and the membership status of our two remaining candidates. 

3 SPECTROSCOPIC FOLLOW-UP OF SDSS 

J121845.69+264831.7 AND SDSS J121856.17+254557.1 

3.1 The WHT observations 

We have obtained optical spectroscopy of both SDSS 
J121845. 69+264831.7 and SDSS J121856.17+254557.1 us- 
ing the William Herschel Telescope (WHT) and the double-armed 
Intermediate dispersion Spectrograph and Imaging System 
(ISIS). These observations were conducted on 2006/02/02 and 
2008/07/25, when excellent conditions prevailed. On both nights 
the sky was very clear with seeing ^0.6-0.9", while additionally 
on 2008/07/25, humidity levels were very low (~10%). For the 
observation of SDSS J121845.69-F264831.7, ISIS was configured 
with a 1.0" slit and data covering ARi3600-5500A were acquired 
using the R300B grating on the blue arm only (A/5A«1200). 
For the observation of SDSS J121856.17+254557.1, the spectro- 
graph was configured with a 0.6" slit and data covering the two 
wavelength ranges, ARi3600-5500A, 6200-7000A, were obtained 
simultaneously using the R300B (A/5ARi2000) and R1200R 
(A/(5A~ 10000) gratings on the blue and red arms respectively. We 
performed a series of short integrations to obtain total exposure 
times of 20 and 45 minutes on SDSS J121845. 69-1-264831.7 and 
SDSS 1121856.17-^254557.1 respectively. The CCD frames were 
debiased and flat fielded using the IRAF procedure CCDPROC. 
Cosmic ray hits were removed using the routine LACOS SPEC 
(van Dokkum 2001). Subsequently, the spectra were extracted 
using the APEXTRACT package and wavelength calibrated by 
comparison with the CuAr-l-CuNe arc spectra taken immediately 
before and/or after the science exposures. The removal of re- 
maining intrument signature from the science spectra obtained on 
the first night was undertaken using an observed and a synthetic 
spectrum (Bohlin 2000) of G191-B2B. Data obtained on the 
second night were corrected using an observation of the bright DC 
white dwarf WD 19 18-1-386. The reduced spectra clearly show both 
objects to be hydrogen rich white dwarfs, ie. WD1216-I-270 and 
WD 121 6-1-260. 



3.2 The model atmosphere calculations 

We have used recent versions of the plane-parallel, hydrostatic, 
non-local thermodynamic equilibrium (non-LTE) atmosphere and 
spectral synthesis codes TLUSTY (v200; Hubeny 1988, Hubeny 
& Lanz 1995) and SYNSPEC (v48; Hubeny, I. and Lanz, T. 
2001, http://nova.astro.umd.edu/) to generate a grid of pure-H syn- 
thetic spectra covering the effective temperature and surface grav- 
ity ranges Tcfi =15000-30000K and log 3=7.5-8.75 respectively. We 
have employed a model H atom incorporating the 8 lowest energy 
levels and one superlevel extending from n=9 to n=80, where the 
dissolution of the high lying levels was treated by means of the 
occupation probability formalism of Hummer & Mihalas (1988), 
generalised to the non-LTE situation by Hubeny, Hummer & Lanz 
(1994). The calculations assumed radiative equilibrium (since our 
test calculations and those of the referee, Detlev Koester, indicate 
that in this effective temperature regime convection transports 1 % 
of the flux and its neglect has a neglible impact on the emergent 
spectrum) and included the bound-free and free-free opacities of 
the H~ ion and incorporated a full treatment for the blanketing ef- 
fects of HI lines and the Lyman —a, — /3 and —7 satellite opac- 
ities as computed by N. Allard (e.g. Allard et al. 2004). During 
the calculation of the model structure the hydrogen line broadening 
was addressed in the following manner: the broadening by heavy 
perturbers (protons and hydrogen atoms) and electrons was treated 
using Allard's data (including the quasi-molecular opacity) and an 
approximate Stark profile (Hubeny, Hummer & Lanz 1994) respec- 
tively. In the spectral synthesis step detailed profiles for the B aimer 
lines were calculated from the Stark broadening tables of Lemke 
(1997). 



3.3 Determination of the effective temperature and surface 
gravity 

As in our previous work (e.g. Dobbie et al. 2006a), comparison 
between the models and the low resolution data is undertaken us- 
ing the spectral fitting program XSPEC (Shafer et al. 1991). In our 
analysis all lines from H-/3 to H-8 are included in the fitting process. 
XSPEC works by folding a model through the instrument response 
before comparing the result to the data by means of a x^ — statistic. 
The best fit model representation of the data is found by increment- 
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Table 2. Alternative designations, effective temperatures, surface gravities, predicted absolute g magnitudes and distance moduUi for SDSS 
J121845.69+264831.7 and SDSS J121856.17+254557.1. 



Object 


Other IDs 




log 9* 




(m-M)j, 


SDSS J121845.69+264831.7 
SDSS J121856.17+254557.1 


LB4, WD1216+270 

Ton 607, CSI+26-12165, WD1216+260 


258271^^^ 
15739lJ«^ 


7Q1+0.03 

8.46l°;°3 


9 94+0.11 

"•"^-0.12 

11.84l°;;^2 


6.531"- 
4-68l°;1^2 



* Formal fit errors. 



1 






SDSS 


J121845. 69+264831. 7 




1 1 






SDSS 


J121856. 17+254557.1 




4 
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m-M 

Figure 4. The estimated distance modulii of WD1216+270 (SDSS 
J121845.69+264831.7) and WD1216+260 (SDSS J121856.17+254557.1) 
based on the observed SDSS photometry and the synthetic photometry from 
Holberg & Bergeron (2006). The distance of the centre (d=85pc) and the 
tidal limits (rtidai=6-4pc) of the cluster are overplotted (solid and dotted 
vertical lines respectively). WD1216+260 resides within the tidal bounds 
of Melotte 111. 

ing free grid parameters in small steps, linearly interpolating be- 
tween points in the grid, until the value of is minimised. Errors 
in the TcffS and log g s are calculated by stepping the parameter 
in question away from its optimimi value and re-determining min- 
imum until the difference between this and the true minimum 
corresponds to lu for a given number of free model parameters 
(e.g. Lampton et al. 1976). The results of our fitting procedure are 
given in Table 2 and shown overplotted on the data in Figure 3. It 
should be noted that the parameter errors quoted here are formal Icr 
fit errors and undoubtedly underestimate the true uncertainties. For 
our subsequent calculations we adopt more realistic levels of un- 
certainty of 2.3% and 0.07dex in effective temperature and surface 
gravity respectively (e.g. Napiwotzki, Green & Saffer 1999). 

3.4 Distance, radial velocity and cluster membership status 

Using the effective temperature and the surface gravity estimates 
listed in Table 2, and the model white dwarf photometry of Berg- 
eron et al. (1995), which has been updated by Holberg & Berg- 
eron (2006), we have derived the absolute g magnitudes (Mg) of 
WD1216-I-270 (SDSS J121845.69-I-264831.7) and WD1216-I-260 
(SDSS J121856.17-I-254557.1). Subsequently, we have determined 
the distance modulii of these two objects, neglecting extinction 
which is minimal along this line of sight. These are plotted in Fig- 
ure 4 (solid bars), along with reasonable estimates for the distance 
and the tidal bounds of Melotte 111 (d=85pc, rtidai=6.4pc). Inspec- 
tion of this plot reveals that while WD1216-I-270 is probably a field 
star residing some way behind the cluster, the distance modulus 
of WD 1216+260 argues rather strongly that it is associated with 
Melotte 111. 

To further scrutinise the membership status of this one remain- 
ing candidate (WD 1216+260), we have acquired a radial velocity 
measurement. We have removed the effects of telluric water vapour 
from the red arm ISIS spectroscopic data using a template absorp- 
tion spectrum. Custom written DDL routines have been used to nor- 
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Figure 5. The results of our fitting of a synthetic profile (thin black fine) 
to the central portions of the observed H-a Balmer Une of WD1216+260 
(thick grey fines). The fluxx units are arbitrary. 



malise both the observed H-a line and the profile from a synthetic 
spectrum corresponding to rcfi=15750K and log 3=8.45. Subse- 
quently, the model has then been compared to the data, where the 
Levenberg-Marquardt algorithm has been used to minimise a x^ 
fit statistic. The result of this process is displayed in Figure 5. The 
line velocity shift, after correction to the heliocentric rest frame, is 
shown in Table 3. The uncertainty in this measurement has been 
estimated using the bootstrapping method of statistical resampling 
(Efron 1982). The gravitational redshift component of the velocity 
shift was then determined using Equation 1, where M and R are 
the mass and radius of the white dwarf in solar units respectively 
and V is the gravitational redshift in kms~ ^ . 

V = 0.635M/i? (1) 

The mass and radius have been determined using modern evo- 
lutionary tracks supplied by the Montreal group (e.g. Fontaine, 
Brassard & Bergeron 2001), more specifically the calculations 
which assume a mixed CO core and thick H surface layer structure. 
The radial velocity is represented by the difference between the 
measured shift of the line and the calculated gravitational redshift. 
It can be seen from Table 3 that the radial velocity of WD1216+260 
is very close to the value expected for a member of Melotte 111 
and entirely consistent within the measurement uncertainties. On 
the basis of the projected location, the proper motion, the distance 
and the radial velocity of WD 121 6+260, we conclude that it is a 
white dwarf member of Melotte 111. 
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Table 3. Additional physical details of WD1216+260 (SDSS J121856.17+254557.1), the white dwarf probable member of Melotte 111. The mass 
and the cooling time have been estimated using the mixed CO core composition "thick H-layer" evolutionary calculations of the Montreal Group (e.g. 

Fontaine. Brassard & Bergeron 2001). 



Object 


MwD 
(M0) 


RwD 
(R0X1O-3) 


rVMellll 


H-a shift GR 
(kms-l) 


rvwD 


Tc 

(Myrs) 


Minit 

(M0) 


SDSS J121856.17+254557.1 


0.902 ± 0.045 


9.26±0.84 


-0.4±0.4 


57.lt°j 61.9±6.4 


-4 8+''-'* 


363l« 


A 77+5.37 
'-0.97 




3.60 3.80 4.00 4.20 4.40 

log,o[Wavelength (A)] 

Figure 6. The observed optical (SDSS u,g,r,i and z) to near-IR (2MASS 
PSC J and UKIDSS GCS K) spectral energy distribution of WD1216+260 
(open squares and error bars). The model photometry of Bergeron et al. 
(1995) as updated by Holberg & Bergeron (2006) is overplotted (solid cir- 
cles), where we have assumed T(,fj=15739K and log g = 8.46 and a dis- 
tance modulus of (m-M)g=4.68. A TLUSTY/SYNSPEC synthetic spectrum 
for Toff = 15739K and log g = 8.46 and scaled to match the observed op- 
tical photometry, is also overplotted (sohd Une). The estimated J and K 
flux levels resulting from the presense of a putative spatially unresolved 
M=O.O36M0 brown dwarf companion are shown (asterisks). 



4 A WHITE DWARF MEMBER OF MELOTTE 111 

4.1 The near-IR energy distribution of WD1216-I-260 

We note that due to its proximity and intrinsic luminosity 
WD1216-I-260 was detected in the 2MASS All Sky Survey (Skrut- 
skie et al. 2006), albeit weakly and only in the shortest wave- 
band, J. It is listed in the Point Source Catalogue as having 
J2MAss=16.79±0. 17. Furthermore, it was detected at K in the 
course of the Galactic Clusters Survey (GCS) component of the 
United Kingdom Infrared Digital Sky Survey (UKIDSS; Lawrence 
et al. 2007). Data extracted from the Wide Field Camera Science 
Archive (WSA; Hambly et al. 2008) indicates that for this object, 
/s:wFCAM=17.08±0.09. UKIDSS J and H band photometry are 
not yet available for the Melotte 111 region of sky. 

We have used data listed in Holberg & Bergeron (2006; SDSS 
ugriz), Cohen et al (2003; 2MASS J) and Hewett et al (2006; 
WFCAM K) to convert the optical and near-IR magnitudes to 
fluxes and have plotted these in Figure 6 (open squares). Model 
DA white dwarf photometry, from Bergeron et al. (1995) as up- 
dated by Holberg & Bergeron (2006), is overplotted (solid circles), 
where we have assumed Tcff=15739K and log g = 8.46 and a dis- 
tance modulus of (m-M)g=4.68. Note that the model K magnitude 
has been moved from the 2MASS system to the Wide Field Camera 
(WFCAM) system using transform equations from Hodgkin et al. 



(2009), here extrapolated slightly to the blue. Additionally we have 
overplotted a TLUSTY/SYNSPEC synthetic spectrum of the ap- 
propriate effective temperature and surface gravity, which has been 
scaled to match the observed optical fluxes (solid line). 

Approximately 25% of DAs in the effective temperature 
regime occupied by WD1216-I-260 are metal-rich (ie. DAZs; Zuck- 
erman et al. 2003) and at least 14% of these are observed, through 
the existence of an infrared excess (e.g. KiUc et al. 2005, Beck- 
lin et al. 2005), to harbour a dust disk which may be associated 
with a former planetary system (Kilic & Redfield 2007; DAZd). 
The discovery of a DAZd in an open cluster, where the population 
age and progenitor mass can be relatively well constrained, would 
be of particular interest. However, the good agreement between the 
observed and the theoretical infrared fluxes means that the current 
data provide no evidence for the existence of a dust disk around 
WD1216-I-260. Moreover, these data do not indicate the presense 
of a close (spatially unresolved) cool low mass companion (e.g. 
Farihi & Christopher 2004, Dobbie et al. 2005, Maxted et al. 2006, 
Burleigh et al. 2006). 

Since WD1216-I-260 is a comparatively massive object with a 
likely total age of T~500Myrs, we are able to set a fairly tight upper 
limit on the mass of a putative cool companion. To do this we first 
moved the mean CIT magnitudes for L and T dwarfs listed in Table 
9 of Vrba et al. (2004) onto the 2MASS (J) and Mauna Kea Ob- 
servatories {K; MKO) sytems using published fits relating spectral 
type and magnitude offsets between photometric systems (Stephens 
& Leggett 2004). We assumed here that i^MKO~J^WFCAM, since 
WFCAM employs the MKO filter set (e.g. Hewett et al. 2006) and 
Leggett et al. (2006) find no compelling evidence for a colour term 
at K between the UFTI MKO system (e.g. Hawarden et al. 2001) 
and the WFCAM system (although their data is rather limited). We 
then constructed the relations, spectral type versus luminosity, ef- 
fective temperature and near-IR magnitudes using the Vrba et al. 
(2004) data, where the COND03 models of Baraffe et al. (2003) 
were used to estimate the upwards offset in luminosity between this 
field dwarf dominated sample and the younger open cluster popu- 
lation. As we were aiming to set an upper limit on the mass of a 
putative companion, we assumed T=600Myrs for Melotte ill (the 
upper bound of the age estimate) and further offset the predicted 
magnitudes by -1-0.4 to account for the dispersion observed in the 
near-IR photometry of L/T dwarfs of a given spectral type (e.g. 
Tinney et al. 2003). Subsequently, beginning with T8 we combined 
progressively earUer and earlier spectral types with the model white 
dwarf photometry until the predicted flux at K matched the 3a up- 
per limit on the observed flux (ifvvFCAM=16.81). This occured at 
spectral type T2.5 (Tefj ~1400K), which for our adopted age corre- 
sponds to M=O.O36M0 . 

Additionally, we inspected a I'xl' (representing a projected 
area of 5000AUx5000AU) WFCAM K band thumbnail image 
centered on WD1216-I-260. We found that this white dwarf is the 
only significant detection within this region. Photometry for other 
catalogued point sources within 30' indicate that a ^3cr detec- 
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tion corresponds to A'^18.65. At the distance of Melotte 111 and 
r=600Myrs this is consistent with M^0.034Mo. Thus we con- 
clude that WD 12 16+260 has no spatially resolved cool companion 
with M<;O.O34M0 within a~2500AU and no spatially unresolved 
cool companion with M<;O.O36M0 . Farihi, Becklin & Zuckermann 
(2005) estimated that ^0.5% of white dwarfs have brown dwarf 
companions with M^O.O2-O.O5M0 within a^5000AU. Thus, this 
result is not entirely unexpected. 

4.2 In the context of the IFMR 

As a member of an open star cluster and having probably evolved 
essentially as a single star, WD 12 16+260 is potentially useful for 
constraining the form of the IFMR. Therefore, we have used the 
Montreal group evolutionary tracks to determine a cooling time, 
where cubic splines have been used to interpolate between points in 
this grid. The results of this process are shown in Table 3. The life- 
time of the progenitor star has been derived by subtracting this cool- 
ing time from the adopted cluster age, where, for the reasons out- 
lined in the Section 1, it has been assumed that r=500±100Myrs. 
Subsequently, we have used cubic splines to interpolate between 
the lifetimes calculated for stars of solar composition by Girardi et 
al. (2000) and have constrained the progenitor mass to the value 
shown in the final column of Table 3. The error we quote in the 
progenitor mass takes into account the spread in recent cluster age 
determinations and the uncertainty in the cooling time of the white 
dwarf, but should be regarded as an illustrative rather than a robust 
estimate. 

We find the location of WD1216+260 in initial mass-final 
mass space to be consistent with the general trend outlined by 
the bulk of white dwarf members of near solar metalicity (ie. 
-0.15^[Fe/H]^+0.15) open clusters (e.g. Kalirai et al. 2008 - 
NGC6819 and NGC7789, Claver et al. 2001 - Hyades and Prae- 
sepe, Casewell et al. 2009 - Praesepe, Kalirai et al. 2005 - 
NGC2099, Williams & Bolte 2007 - NGC6633, Rubin et al. 2008 - 
NGC1039, Dobbie et al. 2009 - NGC2287 and NGC3532, Koester 
& Reimers 1996 - NGC2516 and Dobbie et al. 2006a,b - Pleiades). 
While the mass range for the progenitor star of WD1216+260 is 
broad, primarily due to the substantial uncertainty on the age of 
Melotte 111 (e.g. Salaris et al. 2008), the most probable initial and 
final masses are very similar to those of the two white dwarf mem- 
bers of NGC2287 and the most massive degenerate in NGC1039. 
Therefore, WD 12 16+260 lends some support, albeit rather weak 
given the large uncertainties, to our conclusion in Dobbie et al. 
(2009) that the cluster data appear indicate that the IFMR is less 
steep at Minit >4M0 than in the range 3M0<Minit^4M0. As 
discussed at length in Dobbie et al. (2009), this decrease in the gra- 
dient of the IFMR is expected on theoretical grounds. The form of 
the initial mass versus the core mass at the time of the first ther- 
mal pulse relation (Karakas, Lattanzio & Pols 2002) and a recent 
theoretical IFMR (Marigo & Girardi 2007), reasonably reproduce 
that of the data at Minit>3M0, although admitedly with offsets in 
final mass. Additional and/or better data in the initial mass regime 
Minit <;5M0 would help to consolidate this conclusion. 



5 SUMMARY 

We have utilised SDSS photometry, SuperCOSMOS proper mo- 
tions and optical spectroscopy to search for white dwarf members 
within the central regions of the open cluster Melotte 111. We 
have identified one probable degenerate member, WD1216+260, 



with a mass of Mwd=0.90 ± O.O4M0. The optical to near-IR en- 
ergy distribution of this object, as described by SDSS, 2MASS and 
UKIDSS GCS data, is entirely consistent with that of an isolated 
DA white dwarf. We have set tight upper limits on the mass of a 
putative low mass companion, M<;O.O36M0 (spatially unresolved) 
and M^O.O34M0, (spatially resolved and a^2500AU). Moreover, 
for an adopted cluster age of r=500±100Myrs, we have infered the 
mass of its progenitor star to be Minit=4.77io g7M0. The location 
of WD1216+260 in initial mass-final mass space is consistent with 
a near-monotonic positive correlation outlined by the bulk of open 
cluster white dwarfs and our earlier conclusion that the these ob- 
jects appear to indicate that the IFMR is less steep at Minit>4M0 
than in the range 3M0 ;$Minit ^4M0 , at least for metalicity close 
to the solar value. 
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